Could MACHOS be Primordial Black Holes formed during the QCD Epoch? 
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Observations by the MACHO collaboration indicate that a significant fraction of the galactic halo 
dark matter may be in form of compact objects with masses M ~ 0.5M©. Identification of these 
objects as red or white dwarfs is problematic due to stringent observational upper limits on such 
dwarf populations. Primordial black hole (PBH) formation from pre-existing density fluctuations 
is facilitated during the cosmic QCD transition due to a significant decrease in pressure forces. 
For generic initial density perturbation spectra this implies that essentially all PBHs may form 
with masses close to the QCD-horizon scale, M^ CD ~ 1M©. It is possible that such QCD PBHs 
contribute significantly to the closure density today. I discuss the status of theoretical predictions for 
the properties of QCD PBH dark matter. Observational signatures of and constraints on a cosmic 
solar mass PBH population are also discussed. 



I. PBH FORMATION DURING THE QCD 
EPOCH 



It is long known that only moderate deviations from 
homogeneity in the early universe may lead to abundant 
production of PBH's from radiation [l]]. For a radiation 
equation of state (i.e. p = p/3, where p is pressure and p 
is energy density) there is approximate equality between 
the cosmic Jeans-, Mj, and horizon-, Mh, masses. The 
ultimate fate of an initially super-horizon density fluc- 
tuation, upon horizon crossing, is therefore determined 
by a competition between dispersing pressure forces and 
the fluctuation's self-gravity. For fluctuation overden- 
sities exceeding a critical threshold at horizon crossing 
(Sp/p) hc > S^ D w 0.7 H formation of a PBH with mass 
M pbh ~ Mh results. 

The universe must have passed through a color- 
confinement quantum chromodynamics (QCD) transi- 
tion at cosmic temperature T w lOOMeV. Recent lat- 
tice gauge simulations indicate that the transition be- 
tween a high-temperature quark-gluon phase and a low- 
temperature hadron phase may be of first order || , even 
though such simulations are still plagued by limited res- 
olution and problems to account for finite strange quark 
mass. A first order phase transition is characterized by 
coexistence of high- and low-temperature phase at coexis- 
tence temperature T c . Both phases may exist in pressure 
equilibrium, p^P — p^ but with different and constant (at 
T c ) energy densities, p q ® — p 1 ^ = L, where L is the latent 
heat. During phase coexistence adiabatic expansion of 
the universe causes a continuous growth of the volume 
fraction occupied by hadron phase (1 — f qg ), on the ex- 
pense of quark-gluon phase, such that through the release 
of latent heat the universe is kept at T c . The transition 
is completed when all space is occupied by hadron phase. 

Consider a volume element of mixed quark-gluon- and 
hadron phase- during phase coexistence. Provided a typi- 
cal length scale of the volume element is much larger than 
the mean separation between quark-gluon and hadron 
phase (i.e. the mean hadron- or quark-gluon- bubble 
separation, l s ) one may regard the volume element as ap- 
proximately homogeneous. The average energy density of 



the volume element is (p) — p^ + f qg L and continuously 
varies with the change of f qg , whereas pressure remains 
constant, p = pjf. Upon adiabatic compression of mixed 
quark-gluon/hadron phase there is therefore no pressure 
response, v e Jf — y/ (dp/d(p)) s =0 0. Of course, the 
pressure response may only vanish if thermodynamic 
equilibrium is maintained. For rapid compression time 
scales or small compression amplitudes this may not be 
the case, whereas it is anticipated that approximate ther- 
modynamic equilibrium applies over a Hubble time and 
order unity compression factors. During phase coexis- 
tence the universe is effectively unstable to gravitational 
collapse for all scales exceeding l s . Note that a vanishing 
of v%ff, which was independently discovered by ||, may 
also have interesting non-gravitational effects on density 
perturbations. 

These considerations have led me Q to propose PBHs 
formed during the QCD epoch from pre-existing initially 
superhorizon density fluctuations, such as leftover from 
an early inflationary period of the universe, as a candi- 
date for non-baryonic dark matter. Fluctuations cross- 
ing into the horizon during the QCD epoch experience a 
significant reduction of pressure forces over that regime 
of the fluctuation which exists in mixed phase. Since 
the PBH formation process is a competition between 
self-gravity and pressure forces, and v s — 1/V3 is con- 
stant during most other radiation dominated epochs, the 
threshold for PBH formation should be smaller during 
the QCD epoch than during other early eras, <5^ CD < 
<5^ D . Only a slight favor for PBH formation during 
the QCD epoch may effectively lead to the production 
of PBH on only approximately the QCD horizon mass 
scale M^ CD « 2M Q (T c /100MeV)- 2 . This holds true for 
strongly declining probability distribution functions for 
the pre-existing fluctuation overdensities. For example, 
assuming Gaussian statistics, PBH formation is domi- 
nated for Sp/ p in the range S c and S c + a 2 / 'S c , where a is 
the variance of the Gaussian distribution. This range is 
very small, ij 2 /5 c < 10~ 2 , if PBH mass density is not to 
exceed the present closure density, fl p bh ^ 1- For Q p bh ~ 1 
PBH formation during the QCD epoch is also a very rare 
event with only a fraction ~ 10 -8 of horizon volumes col- 
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lapsing to black holes. 

The possible production of PBH during the QCD 
epoch is not a completely new suggestion. In fact, in 
the mid seventies it was believed that a QCD era was 
characterized by an ever-increasing production of massive 
hadronic resonances. Such a " soft " (i.e. almost pressure- 
less) Hagedorn era was argued to be suspect since over- 
production of primordial black holes seemed likely || . In 
the eighties it was argued that the long-range color force 
could lead to the generation of subhorizon density fluc- 
tuations which in turn could collapse to planetary sized 
PBHs Nevertheless, the simple properties of mixed 
phase during a cosmic first-order transition and their pos- 
sible implications for PBH formation on the QCD horizon 
mass scale have so far been overlooked. 

Currently there are two groups attempting to simulate 
the PBH formation process during a QCD transition with 
aid of a general- relativistic hydrodynamics code ^[j). 
Preliminary results by || verify the reduction of PBH 
formation threshold during the QCD epoch. Assuming 
a bag equation of state and phase transition parameter, 
Lj p^ = 2, we have found a PBH formation threshold re- 
duction, 5® /5f D « 0.77, for fluctuations entering the 
horizon approximately during the middle of the phase 
transition. Note that a canonical bag model with to- 
tal statistical weights of g qg = 51.25 and = 17.25 
for quark-gluon- and hadron- phases [Q, respectively, 
predicts even larger L/ = 2.63, whereas lattice simula- 
tions may favor smaller L |nj . For fluctuations entering 
the horizon during the QCD epoch one typically finds 
the evolution of the fluctuation into two different spatial 
regimes. An inner part of the fluctuation exists in pure 
guark-gluon phase p > p q ® surrounded by an outer part 
existing in pure hadron-phase p < p^. The enhanced 
density in the inner part of the fluctuation assists the 
collapse to a PBH. This is in contrast to PBH forma- 
tion during simple radiation dominated eras, where the 
fluctuation's density distribution is continuous. It has 
also been attempted to derive approximate analytic es- 
timates of the threshold reduction and PBH masses for 
PBH formation during the QCD epoch |l2[ . The model 
predicts that S® CB is minimized for fluctuations entering 
the horizon well before the transition resulting into PBH 
masses considerably smaller than the QCD horizon. 



II. THEORETICAL PREDICTIONS FOR QCD 
PBH SCENARIOS 



It is valuable to advance the initial suggestion of possi- 
ble abundant production of PBH during the QCD epoch 
to a complete and predictive scenario. I outline here to 
which degree this may be accomplished and briefly de- 
scribe the theoretical issues in QCD PBH scenarios. 

Threshold reduction: The bias for forming PBHs 
almost exclusively on the QCD scale is dependent on 
the PBH formation threshold reduction, <5? CD < (5? D . 



Within the context of a bag model equation of state for a 
first order transition preliminary results of numerical sim- 
ulations confirm the proposed threshold reduction. For 
higher order QCD transitions threshold reduction could 
still occur but would have to be verified by using accu- 
rate p{T), p{T), and v s (T) determined from lattice gauge 
simulations. Due to the duration of the QCD transition, 
L/Pc ~ 1; threshold reduction will be of order unity. A 
very accurate determination of <5^ CD is only necessary if 
PBH formation is efficient over a range in dp/p which 
is also of order unity. This is not the case for Gaussian 
statistics of the pre-existing density fluctuations but may 
apply for non-Gaussian statistics. 

Mass function: A crucial prediction of a QCD PBH 
dark matter scenario is the average QCD PBH mass. 
There is seemingly rough agreement of the QCD hori- 
zon mass scale Af^ CD w 2M (T c /lOOMeV)~ 2 and the 
inferred masses of compact objects in the galactic halo 
by the MACHO collaboration, M ~ O.5M . Neverthe- 
less, currently there are large uncertainties in the pre- 
diction for average QCD PBH mass, (M^^). Even in- 
correctly assuming M^ CD = (A/^^ D ), there is a factor 

eight uncertainty in M^ hi ^ depending on if the horizon 
length is taken as radius or diameter of a spherical hori- 
zon volume. The QCD equation of state, order of the 
transition, and transition temperature are as yet not pre- 
cisely determined. The transition temperature may fall 
somewhere in the range 200 MeV > T c > 50 MeV imply- 
ing a factor sixteen uncertainty in A/^ CD , and proba- 
bly equal uncertainty in (M^^ ). Assuming a first or- 
der transition, (M®™) may also depend on L and the 
equation of states above, and below, the transition point. 
An accurate determination of (M®™) requires detailed 
and reliable lattice gauge simulation data. Approximate 
trends may be obtained by using a bag equation of state. 
A PBH mass function, as well as (M^^ D ), is obtained by 
convolving the distribution function for density contrast 
of the pre-existing density perturbations, Sp/p, with a 
scaling relation associating final PBH mass with density 
contrast jl3). The average PBH mass is thus also de- 
pendent on the statistics of the density perturbations. 
Further, it has been shown that resulting PBH masses 
are dependent on the fluctuation shape H|. These un- 
certainties are particularly difficult to remove since they 
require knowledge about the underlying physics creating 
density perturbations, presumably occurring at a scale 
not accessible to particle accelerators. 

Contribution to O p bh: The contribution of QCD 
PBHs to the closure density at the present epoch is 
dependent on the fraction of space which is overdense 
by more than <5^ CD . COBE normalized, exactly scale- 
invariant [n = 1) Gaussian power spectra, imply neg- 
ligible PBH production. Gaussian blue spectra with 
1.37 < n < 1.42 predict Q, pbh in QCD PBHs in the range 
10~ 5 to 10 3 p2[ . Such spectral indices are consistent with 
cosmic microwave background observations H. Never- 
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theless, blue spectra resulting from inflationary epochs 
have been shown to generically be non-Gaussian, skew- 
negative Density perturbations with an exactly 
scale-invariant, COBE normalized power spectrum, but 
with a non-Gaussian, skew-positive distribution tail, may 
yield £l p bh ~ 1- One argument against QCD PBH dark 
matter is the degree of fine-tuning involved for obtaining 
Qpbh ~ 1- 

Accretion around recombination: It is long known 
that black holes may efficiently accrete after the epoch 
of recombination |17|] . Whereas accretion does not ap- 
preciably change the black hole masses, conversion of ac- 
creted baryon rest mass energy into radiation may pro- 
duce substantial radiation backgrounds. The presently 
observed X-ray and/or UV backgrounds may be incom- 
patible with a population of PBHs with mass M p bk > 
10 4 M Q and Q pbh > 0.1 @, A population of M pbh ~ 
1M PBHs with large Q p bh is consistent with the ob- 
served X-ray and/or UV backgrounds. Accretion of 
baryons on PBH shortly before the epoch of recombi- 
nation may produce distortions in the blackbody of the 
cosmic microwave background radiation. PBHs with 
Mpbh ~ 1-W© would accrete at the Bondi rate, with 
Thomson drag inefficient. Tidal interactions between 
the accreting gas and neighboring PBHs would lead to 
the transfer of angular momentum and the formation of 
disks around the PBH. Preliminary results of an inves- 
tigation of PBH accretion before recombination indicate 
that the resulting blackbody distortions would be below 
the current FIRAS limit. 

PBH formation during other epochs: Efficient 
PBH formation during the QCD era may, in principle, 
imply formation of PBHs during other epochs as well. 
For example, during the e + e~-annihilation there is a de- 
crease in the speed of sound which may result in a bias 
to form PBHs on the approximate horizon scale of this 
era. Further, for power spectra of the underlying density 
distribution characterized by n > 1 QCD PBH formation 
may be accompanied by PBH formation at earlier times 
on mass scales M <C M^ CD . It is important to verify 
that such PBHs do not violate observational constraints 




III. OBSERVATIONAL SIGNATURES OF QCD 
PBH DARK MATTER 

Ultimately, only by observational technique the exis- 
tence of a population of QCD PBH may be established. 
It is therefore important to establish the observational 
signatures of QCD PBH dark matter. Particular empha- 
sis is laid on observations which may be performed in the 
not-to-distant future. 

Galactic halo microlensing searches: The recent 
results of microlensing searches for compact, galactic halo 
dark matter by the MACHO collaboration Jl{| provide 
some motivation for QCD PBH dark matter. Low event 



statistics as well as uncertainties about the halo model 
which is to be adopted result in fairly large ranges for 
the average MACHO mass, 0.1 AT© < M< 1M Q , and halo 
dark matter fraction provided by MACHOs, Jm Z 0-2. 
The error bars may be reduced by increasing the number 
of observed microlensing events and observing towards 
several line-of-sights (e.g. towards the Large and Small 
Magellanic Clouds). Nevertheless, it will not be possible 
by only the observational MACHO project to determine 
an accurate mass function. Only in combination with 
follow-up observations, such as by a space interferometry 
sattclitc, degeneracy between MACHO lens mass, dis- 
tance, and projected velocity may be lifted and a mass 
function may be determined. 

Alternative interpretations of the MACHO re- 
sults: The inferred masses of MACHOs are close to 
those of stars, stellar remnants, or brown dwarfs. The 
most straightforward interpretation of the observations 
by the MACHO collaboration are that baryonic objects 
have been detected. However, one has to resort to fairly 
extreme galactic models in order for a characteristic MA- 
CHO mass of M < 0.1M Q to be consistent with the ob- 
servations and for brown dwarfs to remain a viable inter- 
pretation for the lenses. A significant contribution to the 
halo dark matter by red dwarfs seems ruled out by obser- 
vations of the Hubble deep field ^o). Halo white dwarfs 
with halo dark matter fractions exceeding fm ^ 0.1 seem 
also in conflict with observations of the Hubble deep field, 
even though this constraint is dependent on somewhat 
uncertain white dwarf ages and cooling curves In 
addition, it has been argued that the light which would 
be emitted by the progenitors of abundant halo white 
dwarf populations has not been observed in deep galaxy 
surveys p^| . It has been suggested that the lenses re- 
sponsible for the observed microlensing are not within 
the halo, but within a warped or thick galactic disk. Such 
scenarios may possibly be rejected by microlensing obser- 
vations on more than one line of sight. There are other 
more, or less, radical interpretations of the results of the 
MACHO collaboration. It is important, not only for the 
viability of QCD PBH dark matter, to establish, or rule 
out, these alternative interpretations. 

Quasar microlensing: The optical depth for mi- 
crolensing of distant quasars by a cosmic component of 
compact, solar mass objects with O c *~ 1 is remarkably 
large. In fact, a constraint of fl c < 0.2 for a population 
of compact objects with masses M c ~ 1M Q has been de- 
rived from observations of broad line radiation- to con- 
tinuum radiation- flux ratios of ~ 100 quasars p3[ . This 
limit relies on the assumption that most continuum radi- 
ation is emitted from within a compact < O.lpc region in 
the center of the quasar, whereas the broad line radiation 
emerges from a much more extended region around the 
quasar. The limit is independent of the clustering proper- 
ties of the compact objects. There is as yet no conclusive 
model for quasar variability. It has thus been proposed 
that quasar variability is due to microlensing of an f2 c ~ 1 
component of compact objects with M c ~ 10~ 3 M Q pi] . 
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QCD PBH dark matter may therefore be constrained by 
large, homogeneous samples of quasar observations, such 
as expected to result from the Sloan Digital Sky Survey, 
hopefully accompanied by an improved understanding of 
the physics of quasars. 

Gravitational wave detection from PBH bina- 
ries: It has been shown that a fraction 1CP 2 — 10 _1 of 
QCD PBHs may form in PBH binaries |25|. This values is 
in rough agreement with the fraction of binaries observed 
by the MACHO collaboration. Gravitational waves emit- 
ted during PBH-PBH mergers are above the expected 
detection threshold for the LIGO/VIRGO interferome- 
ters when occurring within a distance of ~ 15Mpc. For 
galactic halos consisting exclusively of QCD PBH dark 
matter with M p bh ~ 0.5Mq this implies that up to a few 
mergers per year may be detected by the next generation 
gravitational wave interferometers fl25|] . It is particularly 
encouraging that the gravitational wave signal is sensitive 
to the masses of PBH within the binary. One may hope- 
fully also distinguish between neutron star and black hole 
binaries. Establishing the existence of black holes with 
masses well below the upper mass limit for neutron stars 
may strongly argue in favor of primordial black holes. 

Galactic disk accretion: Limits may be placed on 
galactic halo PBH number densities by the accretion in- 
duced radiation which may be observed when a halo PBH 
passes through the galactic disk in the solar vicinity [^6) . 
Nevertheless, even the ~ 10 s objects which will be ob- 
served within the Sloan Digital Sky Survey will not pro- 
vide sufficient statistics to establish, or rule out, an all 
QCD PBH halo with masses as small as ~ 1M@. 



IV. CONCLUSION 

QCD PBHs may be an attractive dark matter can- 
didate. I have outlined here to which degree accurate 
predictions for the properties of QCD PBH dark matter 
may be made. Most uncertain is the contribution to f2 
of such objects since it relies on knowledge about the un- 
derlying density perturbations on mass scales not acces- 
sible to cosmic microwave background radiation observa- 
tions. Predicting QCD PBH mass functions beyond the 
approximate equality between MACHO masses and the 
QCD horizon mass may improve with detailed numeri- 
cal simulations of the PBH formation process and future 
results of lattice gauge simulations for the QCD equa- 
tion of state. A combination of observational techniques, 
such as galactic microlensing searches, quasar microlens- 
ing searches, and gravitational wave interferometry may 
point towards the abundant existence of such objects. Ul- 
timately, the unambiguous detection of a black hole well 
below the maximum mass for neutron stars may argue 
strongly for its primordial nature. 
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